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Abstract. We propose that classical random matrix models are 
properly viewed as finite difi^erence schemes for stochastic differen- 
tial operators. Three particular stochastic operators commonly arise, 
each associated with a familiar class of local eigenvalue behavior. The 
stochastic Airy operator displays soft edge behavior, associated with 
the Airy kernel. The stochastic Bessel operator displays hard edge 
behavior, associated with the Bessel kernel. The article concludes 
with suggestions for a stochastic sine operator, which would display 
bulk behavior, associated with the sine kernel. 



1. Introduction 

Through a number of carefuUy chosen, eigenvalue-preserving transfor- 
mations, we show that the most commonly studied random matrix distri- 
butions can be viewed as finite difference schemes for stochastic differential 
operators. Three operators commonly arise — the stochastic Airy, Bessel, 
and sine operators — and these operators are associated with three familiar 
classes of local eigenvalue behavior — soft edge, hard edge, and bulk. 

For an example, consider the Hermite, or Gaussian, family of random 
matrices. Traditionally, a random matrix from this family has been defined 
as a dense Hermitian matrix with Gaussian entries, but we show that such 
a matrix is equivalent, via similarity, translation, and scalar multiplication, 
to a matrix of the form 

1 2 

—A + diag_i(xi, . . .,Xn~i) + —= ■ "noise", 

in which A is the n-by-n second difference matrix, diag_]^(a;i, . . . ,a;„_i) is 
an essentially diagonal matrix of grid points, and the remaining term is 
a random bidiagonal matrix of "pure noise." We claim that this matrix 
encodes a finite difference scheme for 
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which is the inspiration for the stochastic Airy operator. (The "noise" term 
will be made precise later.) 

The idea of interpreting the classical ensembles of random matrix theory 
as finite difference schemes for stochastic differential operators was origi- 
nally presented in July 2003 Q^, and the theory was developed in 
The present article contains several original contributions, including firm 
foundations for the stochastic Airy and Bessel operators. 



The standard technique for studying local eigenvalue behavior of a ran- 
dom matrix distribution involves the following steps. (1) Choose a family 
of n-hy-n random matrices, n = 2, 3, 4, ... , (2) Translate and rescale the 
nth random matrix to focus on a particular region of the spectrum, and 
(3) Let n oo. When this procedure is performed carefully, so that the 
eigenvalues near zero approach limiting distributions as n cxd, the limit- 
ing eigenvalue behavior often falls into one of three classes: soft edge, hard 
edge, or bulk. 

The largest eigenvalues of many random matrix distributions, notably 
the Hermite (i.e., Gaussian) and Laguerre (i.e., Wishart) ensembles, display 
soft edge behavior. The limiting marginal density, as the size of the matrix 
approaches infinity, of a single eigenvalue at the soft edge is associated 
with the Airy kernel. Tracy and Widom derived formulas for these density 
functions in the cases /3 = 1, 2, 4, relating them to solutions of the Painleve 
II differential equation. See Figure ITTlT al. Relevent references include 

00 mini 123 El El 

The smallest eigenvalues of some random matrix distributions, notably 
the Laguerre and Jacobi ensembles, display hard edge behavior. The lim- 
iting marginal density of a single eigenvalue at the hard edge is associated 
with the Bessel kernel. Formulas exist for these density functions as well, 
expressible in terms of solutions to Painleve equations. See Figure ^jjb). 
Relevant references include [S] El 1^ ■ 

The eigenvalues in the middle of the spectra of many random matrix 
distributions display bulk behavior. In this case, the spacing between con- 
secutive eigenvalues is interesting. The spacing distributions are associ- 
ated with the sine kernel, and formulas for the density functions, due to 
Jimbo, Miwa, Mori, Sato, Tracy, and Widom, are related to the Painleve 
V differential equation. See Figure ILlf c). Relevant references include 
[3jl,i2_,,l&. 

This article contends that the most natural setting for soft edge behavior 
is in the eigenvalues of the stochastic Airy operator 



(1.1) 
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Figure 1.1. Soft edge, hard edge, and bulk distributions, 
associated with the Airy, Bessel, and sine kernels, respec- 
tively. 



and the most natural setting for hard edge behavior is in the singular values 
of the stochastic Bessel operator 

(1.2) _2V5i- + ^+ 2 ,^^.^^„^ 

ax yjx VP 

A suggestion for a stochastic sine operator, along the lines of p.lH1.2(l . 
is presented at the end of the article. The correct interpretations of the 
"noise" terms in (|1.1|) and (|1.2|l will be specified later in the article, as will 
boundary conditions; see Definitions 13.21 and K14l The parameter /3 has its 
usual meaning from random matrix theory, but now the cases /3 = 1,2,4 
do not seem special. 

Numerical evidence is presented in Figures 11.21 and 11.31 The first com- 
pares histograms of stochastic Airy eigenvalues to the soft edge distribu- 
tions of Figure iLTT a'). and the second compares histograms of stochastic 
Bessel singular values to the hard edge distributions of Figure ^T^b). The 
computations were based on the Rayleigh-Ritz method. They are explained 
in further detail in Sections 17. II and 17.21 

The stochastic Airy, Bessel, and sine operators were discovered by inter- 
preting the classical ensembles of random matrix theory as finite difference 
schemes. We argue that 

(1) When scaled at the soft edge, the Hermite and Laguerre matrix 
models encode finite difference schemes for the stochastic Airy op- 
erator. 

(2) When scaled at the hard edge, the Laguerre and Jacobi matrix 
models encode finite difference schemes for the stochastic Bessel 
operator. 
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Figure 1.2. Least eigenvalue of the stochastic Airy oper- 
ator. In each plot, the smooth curve is a soft edge (Tracy- 
Widom) density, and the jagged curve is a histogram of the 
least eigenvalue from 10^ random samples of the stochastic 
Airy operator. The small positive bias in each histogram 
results from the use of a Rayleigh-Ritz procedure, which 
overestimates eigenvalues. 




(a) /3 = 4, a = 



(b) /3 = 2, a = 



(c) /? = 1, a = 



Figure 1.3. Least singular value of the stochastic Bessel 
operator. In each plot, the smooth curve is a hard edge 
density, and the jagged curve is a histogram of the least 
singular value from 10^ random samples of the stochastic 
Bessel operator. The small positive bias in each histogram 
results from the use of a Rayleigh-Ritz procedure, which 
overestimates singular values. 
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See Section for an overview. Exactly what is meant by "scaling" will 
be developed later in the article. Typically, scaling involves subtracting a 
multiple of an identity matrix and multiplying by a scalar to focus on a 
particular region of the spectrum, along with a few tricks to decompose the 
matrix into a random part and a nonrandom part. The structured matrix 
models introduced by Dumitriu and Edclman Jj and further developed by 
Killip and Nenciu TU' and Edclman and Sutton T play vital roles. 
The original contributions of this article include the following. 

• The stochastic Airy and Bessel operators are defined. Care is taken 
to ensure that the operators involve ordinary derivatives of well 
behaved functions, avoiding any heavy machinery from functional 
analysis. 

• The smoothness of eigenf unctions and singular functions is inves- 
tigated. In the case of the stochastic Airy operator, the fcth eigen- 
function is of the form fk4>, in which fk is twice differentiable and 
is a once differentiable (specifically C"^/^^) function defined by 
an explicit formula. This predicts structure in the eigenvectors of 
certain rescaled matrix models, which can be seen numerically in 
Figure^31 Figure lOl considers analogous results for the stochastic 
Bessel operator. 

• The interpretation of random matrix models as finite difference 
schemes for stochastic differential operators is developed. This ap- 
proach is demonstrated for the soft edge of Hermite, the soft and 
hard edges of Laguerre, and the hard edge of Jacobi. 

Notable work of others includes [2] and Although the stochas- 

tic Airy operator is not explictly mentioned in the large /3 asymptotics of 
Dumitriu and Edclman it appears to play an important role. The sto- 
chastic operator approach has very recently been given a boost by Ramirez, 
Rider, and Virag jE], who have proved a conjecture contained in PIE] 
relating the eigenvalues of the stochastic Airy operator to soft edge behav- 
ior. In addition, they have used the stochastic Airy operator to describe 
the soft edge distributions in terms of a diffusion process. 

The next section reviews necessary background material and introduces 
notation. Sectional provides formal definitions for the stochastic Airy and 
Bessel operators and provides an overview of our results, which are devel- 
oped in later sections. 

2. Background 

Much work in the field of random matrix theory can be divided into two 
classes: global eigenvalue behavior and local eigenvalue behavior. 
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Figure 1.4. A prediction of the stochastic operator ap- 
proach. The entrywise ratio of two eigenvectors of the 
rescaled Hermite matrix model H^^^^ is "smoother" than 
either individual eigenvector. The plots are generated from 
a single random sample of -ff^^ with (3 = 2 and n = 10^. 
A log scale is used for visual appeal. V refers to Matlab's 
gradient function, and indicates two applications of 
the gradient function. See Section 17.31 for details. 



Global eigenvalue behavior refers to the overall density of eigenvalues 
along the real line. For example, a commonly studied distribution on n- 
by-n Hermitian matrices known as the Hermite ensemble typically has a 
high density of eigenvalues near zero, but just a scattering near \/2n by 
comparison. Such a statement does not describe how the eigenvalues are 
arranged with respect to each other in either region, however. 

In contrast, local eigenvalue behavior is observed by "zooming in" on 
a particular region of the spectrum. The statistic of concern may be the 
marginal distribution of a single eigenvalue or the distance between two 
consecutive eigenvalues, for example. Local eigenvalue behavior is deter- 
mined by two factors — the distribution of the random matrix and the region 
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Figure 1.5. Another prediction of the stochastic opera- 
tor approach. The entrywise ratio of two singular vectors 
of the rescaled Jacobi matrix model J^^'^ is "smoother" 
than either individual singular vector. See Section 17.31 for 
details. 



of the spectrum under consideration. For example, the eigenvalues of the 
Hermite ensemble near zero display very different behavior from the eigen- 
values near the edge of the spectrum, at V2n. Conceivably, the eigenvalues 
of a different random matrix may display entirely different behavior. 

Interestingly, though, the eigenvalues of many, many random matrix 
distributions fall into one of three classes of behavior, locally speaking. 
Notably, the eigenvalues of the three classical ensembles of random matrix 
theory — Hermite, Laguerre, and Jacobi — fall into these three classes as the 
size of the matrix approaches infinity. 

In this section, we present background material, covering the three most 
commonly studied random matrix distributions and the three classes of 
local eigenvalue behavior. 

2.1. Random matrix models. There are three classical distributions of 
random matrix theory; Hermite, Laguerre, and Jacobi. The distributions 
are also called ensembles or matrix models. They are defined in this section. 
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Also, joint distributions for Hermite eigenvalues, Laguerre singular values, 
and Jacobi CS values are provided. We use the word spectrum to refer 
to all eigenvalues or singular values or CS values, depending on context. 
Also, note that the language earlier in the article was loose, referring to 
eigenvalues when it would have been more appropriate to say "eigenvalues 
or singular values or CS values." 

2.1.1. Hermite. The Hermite ensembles also go by the name of the Gauss- 
ian ensembles. Traditionally, three flavors have been studied, one for 
real symmetric matrices, one for complex Hermitian matrices, and one 
for quaternion self-dual matrices. In all three cases, the density function is 

(2.1) const X exp tr dA, 

in which (3 = \ for the real symmetric case, [3 — 2 for the complex Hermitian 
case, and /3 = 4 for the quaternion self-dual case. The entries in the upper 
triangular part of such a matrix are independent Gaussians, although the 
diagonal and off-diagonal entries have different variances. 

The eigenvalues of the Hermite ensembles have joint density 

n 

(2.2) const X e-^^^=i^^ J]^ | - Aj/ J]^ dA^ . 

l<i<j<n i—1 

Dumitriu and Edelman extended the Hermite ensembles to all /3 > ^ . 
Below, X ^ Y indicates that X and Y have the same distribution. 

Definition 2.1. The n-by-n (3-Hermite matrix model is the random real 
symmetric matrix 



\/2Gi X(n-l)P 
X(n-l)f) \/2G'2 X(„-2)/3 



in which Gi, . . . , G„ are standard Gaussian random variables, Xr denotes a 
chi-distributed random variable with r degrees of freedom, and all entries 
in the upper triangular part are independent. 

The /3 = 1, 2 cases can be derived by running a tridiagonalization algo- 
rithm on a dense random matrix with density function H2.1(l . a fact first 
observed by Trotter 2A . is the natural extension to general /3, and it 
has the desired eigenvalue distribution. 

Theorem 2.2 (PP). For all fi > 0, the eigenvalues of the (3-Hermite matrix 
model have joint density ijj^.jj)) . 
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As /? — > oo, the /3-Hermite matrix model converges in distribution to 

y/n^ 

y/n-2 
y/n-2 ^^~~3 

^/2 VT 
VT 

This matrix encodes the recurrence relation for Hermite polynomials. In 
fact, the eigenvalues of this matrix are the roots of the nth polynomial, 
and the eigenvectors can be expressed easily in terms of the first ri — 1 
polynomials. See [TOj and ^7j for details. 

2.1.2. Laguerre. The Laguerre ensembles are closely related to Wishart 
matrices from multivariate statistics. Just like the Hermite ensembles, the 
Laguerre ensembles come in three flavors. The /3 — 1 flavor is a distribution 
on real m-by-n matrices. These matrices need not be square, much less 
symmetric. The fi = 2 flavor is for complex matrices, and the /3 = 4 flavor 
is for quaternion matrices. In all cases, the density function is 

(2.3) const X exp ^-^tr dA, 

in which A* denotes the conjugate transpose of A. 

For a Laguerre matrix with m rows and n columns, let a = m — n. It 
is well known that the singular values of this matrix are described by the 
density 

n n 

(2.4) constxe-fSr..A,-Q^^^(a+i)-i -Q |a, _ A^-p J] dA„ 

i—l l<i<j<7i i—1 

in which is the square of the ith singular value. As usual, P — 1 for real 
entries and (3 — 2 for complex entries. 

Dumitriu and Edelman also extended this family of random matrix dis- 
tributions to all /? > and nonintegral a. The notation in this article differs 
from the original notation, instead following jl6| . 

Definition 2.3. The n-by-n j3-Laguerre matrix model, parameterized by 
a > — 1, is the distribution on real symmetric matrices 

X(a+n)f3 X{n-l)l3 

X(a+n-l)l3 X(n-2)P 

X(a+2)P XP 

X(a+l)/3 . 



1 



V3 
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in which Xr denotes a chi-distributed random variable with r degrees of 
freedom, and aU entries are independent. The {n + l)-by-n j3-Laguerre 
matrix models parameterized by a > 0, is 



1 

7^ 



Xnl3 



X(n-l)/3 
X(a+n-2)p X{n~2)0 



X(a+l)/3 X/3 
Xa/3 



with independent entries. 

Notice that {L'^'''){Lf^''')^ and [M^^'^Y {M^-^) are identicaUy distributed 
symmetric tridiagonal matrices. This tridiagonal matrix is actually what 
Dumitriu and Edelman termed the /3-Laguerre matrix model. For more 
information on why we consider two different random bidiagonal matrices, 
see [TC] . 

The [3 — \,2 cases can be derived from dense random matrices following 
the density (|2.3|l . via a bidiagonalization algorithm, a fact first observed by 
Silverstein 15 . Then the general (3 matrix model is obtained by extending 
in the natural way. 

Theorem 2.4 {^). For all (3 > and a > —1, the singular values, 
squared, of L^''^ follow the density 12. 4]) - For all f3 > and a > 0, the 
singular values, squared, of M^''^ follow the density i2.4]) - 

As /3 — > (X), the n-by-n /3-Laguerre matrix model approaches 



(2.5) = 



^/a + n \/n — 1 



a + n — 1 \/n — 2 



VaT2 VT 



%/a+ 1 

in distribution, and the (n + l)-by-n /3-Laguerre matrix model approaches 



(2.6) Af°°'" = 



\/a + n — 1 \Jn—\ 



^ a + n — 2 \fT. 



in distribution. The nonzero singular values, squared, of both of these ma- 
trices are the roots of the nth Laguerre polynomial with parameter a, and 
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the singular vectors are expressible in terms of the first n — 1 polynomials. 
See HSl and U?] for details. 

2.1.3. Jacobi. Our presentation of the Jacobi matrix model is somewhat 
unorthodox. A more detailed exposition can be found in 4 . 

Consider the space of {2n+a+b)-hy-{2n+a+b) real orthogonal matrices. 
A CS decomposition of a matrix X from this distribution can be computed 
by partitioning X into rectangular blocks of size (n + a)-by-n, (n + a)-by- 
{n + a + b), {n + 6)-by-n, and (n + b)-hy-{n + a + b), 



X = 



^11 


Xi2 


X21 


X22 



and computing singular value decompositions for the four blocks. Because 
X is orthogonal, something fortuitous happens: all four blocks have es- 
sentially the same singular values, and there is much sharing of singular 
vectors. In fact, X can be factored as 



(2.7) X = 



C/i 






c 


s ' 




U2 \ 




-s 


c 



in which Ui, U2, Vi, and V2 are orthogonal and C and S are nonneg- 
ative diagonal. This is the CS decomposition, and the diagonal entries 
ci,C2, . . . ,c„ of C are knows as CS values. An analogous decomposition 
exists for complex unitary matrices X, involving unitary Ui, U2, Vi, and 
V2. 

The Jacobi matrix model is defined by placing Haar measure on X. The 
resulting distribution on CS values is most conveniently described in terms 
of Ai = cf , i = 1, . . . , rt, which have joint density 

n n 

(2.8) const X n Af ^"^+'^-'(1 - A,)^(''+^)-^ \{ \K - A,f \{ dK. 

i—l i<j i—1 

The Jacobi matrix model has been extended beyond the real and com- 
plex cases (/3 = 1, 2) to general /3 > 0, first by Killip and Nenciu and later 
by the authors of the present article [U llUlIT^ . The following definition in- 
volves the beta distribution beta(c, d) on the interval (0, 1), whose density 
function is Y^f^^'^^^ ~ xY~^ . 

Definition 2.5. The 2n-by-2n P-Jacobi matrix model J^'°'^, parameter- 
ized by /3 > 0, a > — 1, and 6 > —1, is a distribution on orthogonal matrices 
with a special structure called bidiagonal block form. It is defined in terms 
of random angles 9i, . . . ,9n and (f>i, . . . , <f)n-i from [0, §]. All 2n — l angles 
are independent, and their distributions are defined by 

cos^ 9i ^ heta{^{a + i), §{b + i)) 

cos^ ~ beta(|i, |(a + 6 + 1 + i)). 
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The entries of the /3-Jacobi matrix model are expressed in terms of Ci = 
cos 9i, Si = sin^i, = cos (f>i, and = sin (pi. 







-S2ci 

cis'i 


SnSji—l 


S2S[ 


Si 






-C2C1 

-sis[ 


CnSfi—l 
"Sn— lC„_i 


C2s'i 


Cl . 



Theorem 2.6 f^pi llOj ). Partition the 2n-by-2n f3-Jacobi matrix model into 
four blocks of size n-by-n. The resulting CS values, squared, have density 
function \2.f^) . This is true for all (3 > Q. 

As /3 — > 00, the angles 9i, . . . ,6n and (pi, ... , (pn-i converge in distri- 
bution to deterministic angles 0i, . . . , 0„ and 01, ... , (pn-i, whose cosines 
and sines will be denoted Ci, Si, c[, and s[. It is not difhcult to show 

^» ^ \J a+b+2i ^^'^ ~ \l Becaiisc the angles have deterministic 

limits, the matrix model itself converges in distribution to a fixed matrix 
joo,a,6^ The entries of J~.«.'' encode the recurrence relation for Jacobi 
polynomials. The CS values, squared, of J°°''^'^ are the roots of the nth 
Jacobi polynomial with parameters a, b, and the entries of Ui, U2, Vi, and 
V2 are expressible in terms of the first n — 1 polynomials. See jJHI and ^7] 
for details. 



2.2. Local eigenvalue behavior. The three classes of local behavior in- 
dicated in Figure 11.11 are observed by taking n 00 limits of random 
matrices, carefully translating and rescaling along the way to focus on a 
particular region of the spectrum. This section records the constants re- 
quired in some interesting rescalings. 

For references concerning the material below, consult the introduction to 
this article. Note that much of the existing theory, including many results 
concerning the existence of large n limiting distributions and explicit for- 
mulas for those distributions, is restricted to the cases f3 = 1,2,4. Further 
progress in general (3 random matrix theory may be needed before discus- 
sion of general /3 distributions is perfectly well founded. Although these 
technical issues are certainly important in the context of the stochastic op- 
erator approach, the concrete results later in this article do not depend on 
any subtle probabilistic issues, and hence, we dispense with such technical 
issues for the remainder of this section. 
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2.2.1. Soft edge. Soft edge behavior, suggested by Figure [LTT a). can be 
seen in the right (and left) edge of the Hermite spectrum and the right 
edge of the Laguerre spectrum. 

In the case of Hermite, the fcth largest eigenvalue Xn+i~k{H^) displays 
soft edge behavior. Specifically, —y/2n^^^{Xn+i-k{H^) — \/2n) approaches 
a soft edge distribution as n ^ oo. 

In the case of Laguerre, the kth largest singular value an+i-kiL^'"") 
displays soft edge behavior. Specifically, — 2^/^ri^/^(cr„+i_fc(L''''') — 2-y/n) 
approaches a soft edge distribution as n — > oo. Of course, the fcth largest 
singular value of the rectangular model Af^'" displays the same behavior, 
because the nonzero singular values of the two models have the same joint 
distribution. 

2.2.2. Hard edge. Hard edge behavior, suggested by Figure [TTTT b'! . can be 
seen in the left edge of the Laguerre spectrum and the left and right edges 
of the Jacobi spectrum. 

In the case of Laguerre, the fcth smallest singular value ak^L^'"") dis- 
plays hard edge behavior. Specifically, \/2\J1n + a + \a],{L^'°-) approaches 
a hard edge distribution as n — > cx). Of course, the fcth smallest nonzero 
singular value of Af^'"^ displays the same behavior, because the singular 
values of the two matrix models have the same joint density. 

In the case of Jacobi, the fcth smallest CS value displays hard edge 
behavior. Specifically, let Ckk{J^'°'''') be the fcth smallest diagonal entry in 
the matrix C of 12.7|l . when applied to the /3- Jacobi matrix model. Then 
(2n + a + b+ l)cfefc( J^'"'^) approaches a hard edge distribution as n ^ oo. 

2.2.3. Bulk. Bulk behavior is seen in the interior of spectra, as opposed to 
the edges. Suppose that the least and greatest members of the spectrum of 
an n-by-n random matrix are 0{L„) and 0{Rn), respectively, as n ^ oo. 
Then bulk behavior can often be seen in the spacings between consecutive 
eigenvalues near the point (1 — p)Ln + pRn, for any constant p G (0, 1), 
as n — > oo. This is true for the Hermite, Laguerre, and Jacobi ensembles. 
Because this article does not consider bulk spacings in great detail, the 
constants involved in the scalings are omitted. 

2.3. Finite difference schemes. The solution to a differential equation 
can be approximated numerically through a finite difference scheme. This 
procedure works by replacing various differential operators with matrices 
that mimic their behavior. For example, the first derivative operator can 
be discrctized by a matrix whose action amounts to subtracting function 
values at nearby points on the real line, essentially omitting the limit in 
the definition of the derivative. 
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With this in mind, Vm,n is defined to be the m-by-n upper bidiagonal 
matrix with 1 on the superdiagonal and —1 on the main diagonal, 

-1 1 

V = -11 



The subscripts arc omitted when the size of the matrix is clear from context. 
Up to a constant factor, Vm.n encodes a finite difference scheme for the 
first derivative operator when certain boundary conditions are in place. 

The matrix A„ is defined to bo the symmetric tridiagonal matrix with 
2 on the main diagonal and -1 on the superdiagonal and subdiagonal, 

2 -1 
-1 2 -1 



A„ = 



-1 2 -1 
-1 2 



Note that A„ = V„,„+iV^„_|_i. Under certain conditions, A„ discretizes 
the second derivative operator, up to a constant factor. 

A few other matrices prove useful when constructing finite difference 
schemes. fi„ denotes the n-by-n diagonal matrix with —1, 1,-1,1,... along 
the main diagonal. F„ denotes the n-by-n "flip" permutation matrix, with 
ones along the diagonal from top-right to bottom-left. In both cases, the 
subscript is omitted when the size of the matrix is clear. Finally, the "inter- 
polating matrix" Sm,n = —\^m^ m,n^n provcs useful when constructing 
finite difference schemes for which the domain and codomain meshes inter- 
leave. Sm,n is the m-by-n upper bidiagonal matrix in which every entry on 
the main diagonal and superdiagonal equals ^ . Subscripts will be omitted 
where possible. 



3. Results 

This section defines the stochastic Airy and Bessel operators, briefly 
mentions the stochastic sine operator, and states results that are proved in 
later sections. 

3.1. The stochastic differential operators. 

3.1.1. Stochastic Airy operator. 

Definition 3.1. The classical Airy operator is 
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acting on L^{{0, oo)) functions v satisfying the boundary conditions v{0) = 
0, lim2;^oo v{x) = 0. 

An eigenvalue-eigenfunction pair consists of a number A and a function 
V such that A°°v = Xv. The complete eigenvalue decomposition is 

(^°°-Afc)Ai(-x-Afe) = 0, 

for fc = 1,2,3,..., in which Ai denotes the unique solution to Airy's equa- 
tion f"{x) = xf{x) that decays as a; ^ oo, and Afc equals the negation 
of the fcth zero of Ai. As typical with Sturm-Liouville operators, .4°° 
acts naturally on a subspace of Sobolev space and can be extended to all 
i^((0, oo)) functions satisfying the boundary conditions via the eigenvalue 
decomposition. 

Intuitively, the stochastic Airy operator is obtained by adding white 
noise, the formal derivative of Brownian motion, 

(3.1) -^ + - + ;^^'(-)- 

However, white noise sometimes poses technical difhculties. To avoid these 
potential difficulties, we express the stochastic Airy operator in terms of a 
conjugation of a seemingly simpler operator, i.e., by changing variables. 

Definition 3.2. Let /3 > 0, let B{x) be a Brownian path on (0,oo), and 
let 

(3.2) (j){x) = exp B{x)dx 

The stochastic Airy operator acts on functions v{x) = f{x)(f>{x) satis- 
fying the boundary conditions v{0) — 0, limx^oa v{x) = 0. It is defined 
by 

[A^ifm^) = m ■ ([-4°°/](:r) - A^B{x)r{x) - ^^B{xff{x)^ , 
or, to abbreviate, 

An eigenvalue-eigenfunction pair consists of a real number A and a func- 
tion V such that A^v = Xv. 

Note that 4'{x) is defined in terms of a Riemann integral of Brownian 
motion, which is continuous. This is not a stochastic integral, and nothing 
like an Ito or Stratonovich interpretation must be specified. 

When /3 = oo, the stochastic Airy operator equals the classical Airy 
operator. When /3 < oo, equations (|3.1|l and (|3.3|) are formally equivalent. 
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To see this, apply to v = fcj), proceeding formally as follows. Combining 



and 



/9 ^ 



and 



yields 



The stochastic Airy operator acts naturally on any function of the form 
/0, in which / has two derivatives. Also, the Rayleigh quotient defined by 

{v,A^v) _ li^Tdx + J xv^dx + v^dB 
{v,v) J v^dx ' 

is well defined and does not require an Ito or Stratonovich interpretation if 
V is deterministic, decays sufficiently fast, and is sufficiently smooth, say, 
if it has a bounded first derivative. See 

3.1.2. Stochastic Bessel operator. 

Definition 3.3. The classical Bessel operator with type (i) boundary con- 
ditions, parameterized by a > —1, is the operator whose action is 

dx y'X 

acting on functions v satisfying 

(3.4) type (i) b.c.'s: v{l) = and ij;^v){0) ^ 0. 

We will abuse notation and also denote the classical Bessel operator with 
type (a) boundary conditions by J7^. The action of this operator is also 
J'^ — —2^/x-j^ + and it is defined for all a > —1, but its domain 
consists of functions v satisfying 

(3.5) type (ii) b.c.'s: (jr^^)(0) - and (jr^^)(l) = 0. 
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The adjoint of the classical Bessel operator (with either type (i) or type 
(ii) boundary conditions) has action 

ax y/X 

The singular value decompositions are defined in terms of the Bessel func- 
tions of the first kind ja by 

type (i) b.c.'s: J^[ja{crkVx)] = (Tkja+iifJkVx) 

< (Ti < (72 < CT3 < • • • zeros of ja 

and 

type (ii) b.c.'s: JT\x''l'^\ = 

Ja^baicTkV^)] = CTkja+licrkVx) 
iJrTija+licTkVx)] = Crkja{<JkV^) 

< (Ji < (72 < fa < • • • zeros of ja+l- 

The purposes of the boundary conditions are now clear. The condition at 
X = 1 produces a discrete spectrum, and the condition at a; = eliminates 
Bessel functions of the second kind, leaving left singular functions that are 
nonsingular at the origin. 

Intuitively, the stochastic Bessel operator is obtained by adding white 
noise to obtain 

(3.6) - JT + -^B\ 

with adjoint (jf)* = {J^)* + However, the following definition, 

which avoids the language of white noise, offers certain technical advan- 
tages. 

Definition 3.4. Let a > — 1 and /? > 0, let B{x) be a Brownian path on 
(0, 1), and let 



(3.7) ^{x) 



exp(--L£.-V.,5(, 



The stochastic Bessel operator, denoted J7^, has action 
or, to abbreviate, 

(3.8) J^^i^J^r^. 

Either type (i) or type (ii) boundary conditions may be applied. (See (|3.4|l 
and 
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The adjoint is 

The function ij: involves a stochastic integral, but because the inte- 
grand is smooth and not random, it is not necessary to specify an Ito 
or Stratonovich interpretation. 

Note that when /3 = oo, the stochastic Bessel operator equals the clas- 
sical Bessel operator. When (3 < oo, equations 1)3. 6|l and H3.8|l are formally 
equivalent. To see this, apply J'^ tov — f(t>, proceeding formally as follows. 



= -2V^ f \=x-^l'^B'\ V - 2y/^v' + -^v = J^v + -^B'v. 

The stochastic Bessel operator acts naturally on any function of the form 
fijj for which / has one derivative, assuming the boundary conditions are 
satisfied. Its adjoint acts naturally on functions of the form gip^^ for which 
g has one derivative and the boundary conditions are satisfied. 

Sometimes, expressing the stochastic Bessel operator in Liouvillc normal 
form proves to be useful. The classical Bessel operator in Liouvillc normal 
form, denoted J7^, is defined by 

with either type (i) or type (ii) boundary conditions. The singular val- 
ues remain unchanged, while the singular functions undergo a change of 
variables. The SVD's are 

(3.9) type (i) b.c.'s: J^[Vxja{crkx)] = ak\/xja+i{<7kx) 

{J^y [\/xja+l((7kXy\ = (Jk\/xja{(Tkx) 

< fJi < (72 < 0-3 < • • • zeros of ja 

and 

(3.10) type (ii) b.c.'s: j;^[x''+^'^] = 

J!^[Vxia(.'ykXy\ = akVxja+licTkx) 

{J^)*[Vxja+i{crkx)] = crkVxja{o-kx) 
< fJi < (T2 < (73 < • • • zeros of ja+l- 

Note that although the change of variables to Liouvillc normal form af- 
fects asymptotics near 0, the original boundary conditions still serve their 
purposes. 
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Definition 3.5. Let a > —1 and /? > 0, and let B{x) be a Brownian path 
on (0, 1). The stochastic Bessel operator in Liouville normal form, denoted 
J7jf , has action 



or, to abbreviate, 



ip is defined in l|3.7|l . Either type (i) or type (ii) boundary conditions may 
be appHed. 

This operator acts naturally on functions of the form fip^ for which / 



is once differentiable. It is formally equivalent to ^7^°° + y ^^i?: 

3.1.3. Stochastic sine operator. The last section of this article presents 
some ideas concerning a third stochastic differential operator, the stochastic 
sine operator. This operator likely has the form 

__d_ 



"noise" 


"noise" 


"noise" 


"noise" 



dx 

Key to understanding the stochastic Airy and Bessel operators are the 
changes of variables, in terms of cf) and tp, respectively, that replace white 
noise with Brownian motion. No analogous change of variables has yet been 
found for the stochastic sine operator, so most discussion of this operator 
will be left for a future article. 

3.2. Random matrices discretize stochastic differential operators. 

Much of the remainder of the article is devoted to supporting the following 
claims, relating the stochastic differential operators of the previous section 
to the classical ensembles of random matrix theory. The claims involve 
"scaling" random matrix models. This is explained in Sections |S1 and IHl 

3.2.1. Soft edge ^ stochastic Airy operator. 

Claim 3.6. The Hermite matrix model, scaled at the soft edge, encodes a 
finite difference scheme for the stochastic Airy operator. See Theorems l57 
and I 



Claim 3.7. The Laguerre matrix models, scaled at the soft edge, encode 
finite difference schemes for the stochastic Airy operator. See Theorems 
O and El 
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3.2.2. Hard edge <-> stochastic Bessel operator. 

Claim 3.8. The Laguerre matrix models, scaled at the hard edge, encode 
finite difference schemes for the stochastic Bessel operator. L^'^^^ encodes 
type (i) boundary conditions, and M^!^^ encodes type (ii) boundary con- 
ditions. Sec Theorems EH EHliniZI and EH 

Claim 3.9. The Jacobi matrix model, scaled at the hard edge, encodes 
a finite difference scheme for the stochastic Bessel operator in Liouville 
normal form with type (i) boundary conditions. See Theorems 15.101 and 



3.3. Eigenvalues/singular values of stochastic differential opera- 
tors. Based on the claims in the previous section, we propose distributions 
for the eigenvalues of the stochastic Airy operator and the singular values 
of the stochastic Bessel operators. 

Conjecture 3.10. The kth least eigenvalue of the stochastic Airy operator 
follows the kth soft edge distribution, with the same value for (3. 

The conjecture now appears to be a theorem, due to a proof of Ramirez, 
Rider, and Virag |14| . 

Conjecture 3.11. The kth least singular value of the stochastic Bessel op- 
erator with type (i) boundary conditions follows the kth hard edge distribu- 
tion, with the same values for /? and a. With type (ii) boundary conditions, 
the hard edge distribution has parameters (3 and a + 1 . 

The conjecture should be true for both Jj^ and jj^ . 



This section establishes relations between various random matrices which 
will be useful later in the article. The identities are organized according to 
their later application. This section may be skipped on a first reading. 

4.1. Identities needed for the soft edge. 

Lemma 4.1. Let H^^ = [hij] be a matrix from the n-by-n (3-Hermite 
matrix model. (3 may be either finite or infinite. Let D be the diagonal 
matrix whose {i,i) entry is {n/2)^^^^^''^'^YVkJ^if^k,k+i- Then if f3 < oo, 
DH^D~^ has distribution 



4. Some matrix model identities 



1 



y2Gi 





(4.1) 




^/2G 



"n 
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with all entries independent, while if (3 — oo, DH°°D ^ equals 



(4.2) 



1 

71 





n-l 



n-2 



2 

■spn. 




1 



^ 



Lemma 4.2. Lef L^^'"" = (Zij) 6e a matrix from the n-by-n (3-Laguerre 
matrix model. [3 may be finite or infinite. Let D be the 2n-by-2n diagonal 

matrix whose (i,i) entry is '^-^^'^"'"^^^ 111=1^ ^fcfc nfc=i ^''^^^ Ik.k+i, cmd let 
P = {pij) be the 2n-by-2n "perfect shuffle" permutation matrix, 



Ptj = 



1 j = 2i — 1 or i ~ 2{i — n) 
otherwise 



Then, if (3 < oo, 
(4.3) 



DP' 



(L'3.-)^ 



PD'r 



is the 2n-by-2n random matrix with independent entries 



1 

7^ 





/^A(a+n)/3 















an 



d if 13 = oo, the matrix of \4--^ equals the 2n-by-2n matrix 







n-l 





g+n— 1 



n-2 



Lemma 4.3. Let — {niij) be a matrix from the {n + l)-by-n (3- 

Laguerre matrix model. (3 may be finite or infinite. Let D be the diago- 
nal matrix whose {i, i) entry equals n"^*"-'-)/^ Tli—i^ fTT-kk 111=1 ^'^^^ ™fc+i,fc; 
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and let P ~ (Pij) be the {2n + l)-by-{2n + 1) perfect shuffle matrix, 

'l j ^2i-l or j = 2{i- {n + l)) 
otherwise 

Then, if (3 < go, 

Ml^-" 

is the (2n + l)-by-{2n + 1) random matrix with independent entries 



(4.4) 



DP' 



PD~^ 



1 

7^ 



















/^A(a+n-2)/3 



with all entries independent, and if (3 = oo, the matrix of ^4-4)) equals the 
(2n + l)-by-{2n + 1) matrix 





^ V" 

g+n— 1 Q 
rt-1 



a+n-2 



4.2. Identities needed for the hard edge. The remaining identities are 
derived from the foUowing two completely trivial lemmas. The first oper- 
ates on square bidiagonal matrices, and the second operates on rectangular 
bidiagonal matrices. 

Lemma 4.4. Let A — (aij) and B = (bij) be n-by-n upper bidiagonal 
matrices with the same sign pattern and with no zero entries on the main 
diagonal or superdiagonal. Set 

92i-i = -log|a„;| + log|6ri|, I = 1, . . . ,n 

52i = log|a,j,,,+i| - log z= l,...,n- 1, 

and, for i = 1, . . . , 2n, let di — X]fe=i ^ 9k- Then 

with Deven = diag(d2, ^4, ■ ■ • , c?2n) and Dodd = diag(di, da, . . . , d2„_i). 
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Proof. The {i, i) entry of e^<=™"Be~^°i<i equals e''2*6jie"'^2.-i ^ 5..g-32i-i = 
&jz||f7| = an. The + 1) entry equals e'^'"bi^i+ie~'^^''+^ = 6^,^+16^^' = 
bi^i+i — Qi.i+i- All other entries are zero. □ 

Lemma 4.5. Let A = (aij) and B = (bij) be n-by-{n + l) upper bidiagonal 
matrices with the same sign pattern and with no zero entries on the main 
diagonal or superdiagonal. For i — 1, . . . ,n, set 521-1 = — log |aii| +log \bii\ 
and g2t = log |ai^i+i| - log |6j,i+i|, and for i = 1, . . . ,2n + 1, let di = 
J2k=i 9k- Then A = (,d„,^q^-Dm^ yj^}^ Deven = diag(d2, 0^4, • ■ • , d2n) and 
Dodd = diag(di,d3, . . . ,d2«+i). 



Proof. The (i, i) entry of e^^-^Be"^""''' equals e'^^'&iie"''^'-! = h 
.Kil 



= flij. The entry of e^-™"Se~-°°<id equals e^'^'fei^.+ie-'^^.+i 



= bi^i+i =^ All other entries are zero. □ 

The lemmas immediately establish the following three identities. Con- 
sult Section 1^21 for the definitions of and F. 

For the square /3-Laguerre matrix model, we have 



(4-5) Jtz ^ —FQiL^-^fnF 



-^^j^Fn{L°°'''fnF^ e- 



with 



-Devon = diag((i2, ^4, ■ • ■ , <i2n), -Dodd = diag(di, ds, . . . , d2„-i) 

^k—i 



521-1 -log(^X(a+i)/3) + log Va + j 
52i ~ log(^Xj/3) - logVi. 

5i, ■ ■ • 7 52n-i are independent. 

For the rectangular /3-Laguerre matrix model, we have 



2 



y (2n + a + 1) ^ 



(2n + a + 



e^— ( Wt:^ ^—-^FMM^n'^^n+iFn+^ e"^" 
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with 

Devon = diag(d2, d^, . . . , d2«), -Dodd = diag(di, da, . . . , d2«+i) 
= Efe=j 9k 



921-1 l0g(^X(a+j-l)/3) + log Va + I - 1 

g2i ~ log(^x»/3) - log 

S'l) • • • ) 52n are independent. 

Finally, for the bottom-right block B22^'^ of the /3-Jacobi matrix model 
jl3,a,b^ we have 



with 



^&TTFT^'^S2'^'"^' 

V(2n + a + 6+ ^2 J 

Devon = diag((i2, ^4, • ■ • , rf2n), -Dodd = diag(di, rfs, . . . , d2n-l) 

dv-^2n— 1 
i = l^k=i 9k 

gi (logci - logci) 

g2^ ^ (logs, - log Si) + (logC- - logC-) 

ff2i-i (logc, - logc,) - (logs^_i - logs-_i) for i > 2 

cos^ 9i ^ beta(|(a -4- i), + 



COS'^ I 



2 ' 

beta(fi,f(a + 6+l + ?;)). 



- 2 ' 2 " 

The random angles 0i, . . . , 0„ and 0i, . . . , (fin-i are independent, and their 
cosines and sines are denoted by Ci, Si, c[, and s[, as usual. The constants 
Cj, Sj, c^, and s[ are introduced after the definition of the /3-Jacobi matrix 
model in Section r2. 1.31 

5. Zero temperature matrix models as finite difference 

SCHEMES 

As seen in Section l2.1l the Hermite, Laguerre, and Jacobi matrix mod- 
els approach nonrandom limits as /3 — > cx). We call these matrices "zero 
temperature matrix models" because of the well known connection with 
statistical mechanics. 

By appropriately transforming the zero temperature matrix models — via 
operations such as translation, scalar multiplication, similarity transform, 
and factorization — we can interpret them as finite difference schemes for the 
classical Airy and Bessel operators. This approach anticipates analogous 
methods for the P < oo case. In short, /? = oo matrix models discretize 
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nonrandom operators, and (3 < oo matrix models discretize stochastic op- 
erators. 

5.1. Soft edge. 

5.1.1. Hermite ~* Airy. 

Definition 5.1. Let h — n^^^^. The n-hy-n oo-Hermite matrix model 
scaled at the soft edge is 

in which DH^D^^ is the matrix of (|4.2|l . 

Note that "scahng at the soft edge" modifies eigenvalues in a benign 
way. The translation and rescaling are designed so that the smallest k 
eigenvalues of -ff^ft approach distinct limits asn ^ oo. (The largest eigen- 
values of DH°"D''^ are first pulled toward the origin, and then a scalar 
factor is applied to "zoom in." The scalar factor is negative to produce an 
increasing, as opposed to decreasing, sequence of eigenvalues starting near 
zero.) 

The following theorem interprets iJ^fj as a finite difference scheme for 

,2 

the classical Airy operator A°° — + x on the mesh Xi = hi, i = 

1, . . . , n, with mesh size h = n"^^^. 

Theorem 5.2. For all positive integers n, 

^Toft"^ ■^A + diag_i(a;i,a;2,...,a:„_i), 

in which h = n~^/^ , Xi — hi, and diag_]^(a;i, 2:2, . . . ,a;„_i) is the n-hy-n 
matrix with xi,X2, ■ ■ ■ ,Xn^i on the subdiagonal and zeros elsewhere. 

Furthermore, for fixed k, the kth least eigenvalue of FI'^j^ converges to 
the kth least eigenvalue of A°^' as n—^ 00, 

Proof. The expression for H"^^^ is straightforward to derive. For the eigen- 
value result, recall that the fcth greatest eigenvalue of H°° is the fcth right- 
most root of the nth Hermite polynomial, and the fcth least eigenvalue of 
A°" is the fcth zero of Ai, up to sign. The eigenvalue convergence result 
is exactly equation (6.32.5) of (The recentering and rescaling in the 
definition of i?^ft is designed precisely for the purpose of applying that 
equation.) □ 

It is also true that the eigenvectors of i?^ft discretize the cigcnfunctions 
of .4°°. This can be established with well known orthogonal polynomial 
asymptotics, specifically equation (3.3.23) of ^7]. We omit a formal state- 
ment and proof for brevity's sake. 
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5.1.2. Laguerre Airy. 

1/3 

Definition 5.3. Let h = (2n) The 2n-by-2 n oo -Laguerre matrix 

model scaled at the soft edge is 

in which Dl is the matrix D of Lemma |4.2I (with P — oo) and Pl is the 
matrix P of the same lemma. The {2n-\- l)-by-(2n+ 1) oo-Laguerre matrix 
model scaled at the soft edge is 

PmDiI - V2h-^/^?j , 

in which Dm is the matrix D of Lemma 14.31 (with (3 = oo) and Pm is the 
matrix P of the same lemma. 

L^f^ and M^^^ encode finite difference schemes for the classical Airy 

,2 

operator A°° = —-^ + x. 

Theorem 5.4. Make the approximations 

^Toft = + diag_i(a;i,a;2, . . ■,X2n-i) + El 

and 

^To'ft = + diag_i(a;o, xi, a;2, ■ • • ,X2n-i) + Em, 

with h — (2n)^^/^, Xi — hi, and A\a,g_-^{di, . . . .d„i) denoting the square 
matrix with di, . . . , c?„j on the subdiagonal and zeros elsewhere. Then 
and Em are zero away from the subdiagonal, and their entries are uniformly 
Oih) as n oo. In the special case a = El equals the zero matrix, 
and in the special case a — j, Em equals the zero matrix. 

Furthermore, for fixed k, the kth least eigenvalue of L^^: and the kth 
least eigenvalue of M^'j^ converge to the kth least eigenvalue of the classical 
operator as n oo. 

Proof. For odd j, the (j + entry of El equals — /i(2a+ 1), and every 
other entry El equals zero. For even j, the + 1, j) entry of Em equals 
—h{2a — 1), and every other entry of Em equals zero. For the eigenvalue 
result, check that the kth greatest eigenvalue of 14.3(1 . resp., (|4.4() . equals 
the kth greatest singular value of L°°'°- , resp., M°°'°- , and that this value is 
the square root of the kth rightmost root of the nth Laguerre polynomial 



J oc.a 

^soft ^ 



DlPI 











soft 



T ( DmPL 



M°°^'' 

\,roo,a\T Q 
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with parameter a. The eigenvalue convergence then foUows from equation 
(6.32.4) of jl7| . concerning zero asymptotics for Laguerre polynomials. □ 

Also, the fcth eigenvector of i^f^, resp., M^^^, discretizes the fcth eigen- 
function of via (3.3.21) of 

5.2. Hard edge. 

5.2.1. Laguerre — > Bessel. 

Definition 5.5. The n-by-n oo-Laguerre matrix model scaled at the hard 
edge is 



in which h = — r—rr- F and f2 are defined in Section [231 

Note that "scaling at the hard edge" only modifies singular values by a 
constant factor. This factor is chosen so that the k smallest singular values 
approach distinct limits as n ^ oo. 

The next theorem interprets ^J^^d ^ discretization of the classical 
Bessel operator ^7^°° = + ^ with type (i) boundary conditions. 

Domain and codomain vectors are interpreted on interwoven submeshes. 
The combined mesh has size h — 2n+a+i ^'^^ gy:'iA points Xi = h{a + 
i), i = 1, . . . ,2n. Domain vectors are interpreted on the submesh X2i-i, 
i = 1, . . . ,n, and codomain vectors are interpreted on the submesh X2i, 
i — 1^ . . . ,n. 

Theorem 5.6. Let h and Xi he defined as in the previous paragraph, and 
make the approximation 

^Tard = -2diag(v^, . . . , v^)(^V) 

+ adiag f — • • ■ , / ) S + E, 



with S defined as in Section \2. 51 Then the error term E is upper bidiagonal, 
and the entries in rows , . . . , n of E are uniformly 0{h), for any fixed 
e>0. 

Furthermore, for fixed k, the kth least singular value of L^^^^ approaches 
the kth least singular value of with type (i) boundary conditions as 
n oo, 

'^kiL'^ard) with type (i) h.c.'s). 

Proof. The {i,i) entry of E equals j[y/x2i + ha — \\JX2i — §2^2/^^. By a 
Taylor series expansion, \Jx^ha = y/x + ^x~^/'^ + 0{h'^), uniformly for 
any set of x values bounded away from zero. This implies that the (i,z) 
entry of E is 0{h), for any sequence of values for i bounded below by 
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as n ^ cxD. The (i, z + 1) entry of E equals — |- X2i — ha + \ ^/x2i — f 2:2/^^, 
from which similar asymptotics follow. For the singular value result, recall 
that the fcth least singular value of L°"'°- is the square root of the fcth least 
root of the Laguerre polynomial with parameter a and that the /cth least 
singular value of ^7^°° is the fcth positive zero oi ja^ the Bessel function of 
the first kind of order a. The convergence result follows immediately from 
(6.31.6) of ini. □ 

In fact, the singular vectors of ij^j.^ discretize the singular functions of 
with type (i) boundary conditions as well. This can be proved with 
(3.3.20) of ^13. 

Analogous results hold for the rectangular /3-Laguerre matrix model. 

Definition 5.7. The 7i-by-(n+ 1) 00-Laguerre matrix model scaled at the 
hard edge is 

in which h = — r—rr- F and are defined in Section [TS 

The next theorem interprets Af^^^ as a finite difference scheme for the 
classical Bessel operator = + with type (ii) boundary 

conditions. Domain and codomain vectors are interpreted on interwoven 
submeshes. The combined mesh has size h — 2n+a+i ^'^^ gT:'\d points 
Xi = h{a — 1 + i), i — l,...,2n + l. Domain vectors are interpreted on the 
submesh X2i-i, i = l,...,n+l, and codomain vectors are interpreted on 
the submesh X2i, i = 1, . . . ,n. 

Tlieorem 5.8. Let h and Xi he defined as in the previous paragraph, and 
make the approximation 

MZrd = -2 diag(V^, Vil, . . . , Vi^)(^V„,„+i) 



+ (a - 1) diag ( — ^ ) S'„,„+i + E, 

1X2 Jxl JX2^ 



with S defined as in Section \2. ^ Then the error term E is upper bidiagonal, 
and the entries in rows , . . . , n of E are uniformly 0{h), for any fixed 
e > 0. 

Furthermore, for fixed k, the kth least singular value of M^^^^ approaches 
the kth least singular value of with type (ii) boundary conditions, 

f^fe(^Od) <^k{Jo^i With type (ii) h.c.'s). 

Proof. The {i, i) entry of E equals j^^J X2i + h{a — l) — \s/X2i— ^°2^'' ^2/^^- 
By a Taylor series expansion, ^Jx + h{a - 1) = ^/x+ + 0{h'^), 

uniformly for any set of x values bounded away from zero. This implies that 
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the {i, i) entry of E is 0{h), for any sequence of values for i bounded befow 
by [-1] as n ^ oo. The {i, i + 1) entry of E equals —j^\/x2i — h{a — 1) + 
j^'^X2i — ^^'"■^2/^^' from which similar asymptotics follow. For the singular 
value result, the proof of Theorem 15.61 suffices, because ij^^'^ and M^^"^ 
have exactly the same nonzero singular values, as do with type (i) 
boundary conditions and J^^i with type (ii) boundary conditions. □ 

Also, the singular vectors of M^^'^ discretize the singular functions of 
J^^i with type (ii) b.c.'s, although we omit a formal statement of this fact 
here. 



5.2.2. Jacobi — > Bessel. 

Definition 5.9. The n-by-n 00-Jacobi matrix model scaled at the hard edge 
is 

joo.a.b 1 T71 T->oo,a,b T71 
•^hard = J^^^22 ^ ■• 

in which h = — , ^, , , , and _B?S'"'^ is the bottom-right block of the 2n-bv- 

2n+a+o+l o J 

2n 00-Jacobi matrix model J°°^°-^^ , F and Vl are defined in Section 

The next theorem interprets J^^^'^^^ as a discretization of the classical 
Bessel operator in Liouville normal form, = —^ + {o- + \)'^ with type 
(i) boundary conditions. Domain and codomain vectors are interpreted on 
interwoven meshes. The combined mesh has size h = 7^ — , , , -, and grid 

2n+a+o+l ^ 

points Xi ~ h{a + b + i),i~l,..., 2n. Domain vectors are interpreted on 
the mesh X2,;-i, i = 1, . . . , n, and codomain vectors are interpreted on the 
mesh X2i, i — 1, . . . ,n. 

Theorem 5.10. Let h and Xi be defined as in the previous paragraph, and 
make the approximation 

2h V 2/ Va;2 X4 X2n J 

with S defined as in Section \2. 31 Then the error term E is upper bidiagonal, 
and the entries in rows , . . . , n are uniformly 0{h), for any fixed e > 0. 

Furthermore, for fixed k, the kth least singular value of J^^^'f' approaches 
the kth least singular value of J"^ with type (i) boundary conditions as 
n 00, 

'^kiJZrd'') <^kiJ^ With type (i) b.c.'s). 
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Proof. The entry of J^j.'^'' equals 

1 I a+i / a+b+i 
h\l a + b + 2i\j a + b+l + 2{i-l) 

= ^{X2^ + h{a - b)y^^X-l^\x2r + h{a + b)f'\x2, - hy^'^ 

Rewriting this expression as 

it is straightforward to check that the entry is 577 + (a + ^) • 5^ + 0{h), 
uniformly for any sequence of values i such that X2i is bounded away from 
zero as n — > 00. The argument for the superdiagonal terms is similar. 

For the singular value result, note that the CS values of J°°^°-^'> equal the 
singular values of its bottom-right block, and that these values, squared, 
equal the roots of the nth Jacobi polynomial with parameters a, b. Also 
recall that the fcth least singular value of with type (i) boundary con- 
ditions is the fcth positive zero of ja, the Bessel function of the first kind of 
order a. The rescaling in the definition of J^^d'' designed so that equa- 
tion (6.3.15) of TT| may be applied at this point, proving convergence. □ 

It is also true that the singular vectors of Jj^j^j discretize the singular 
functions of with type (i) boundary conditions. 

As presented here, the theorem only considers the bottom-right block 
of J°°'°-'^^ but similar estimates have been derived for the other three 
blocks JHI ■ Briefly, the bottom- right and top- left blocks discretize and 
(J^)*, respectively, while the top-right and bottom- left blocks discretize 
J'j^ and {J{^)* , respectively, all with type (i) boundary conditions. 

6. Random matrix models as finite difference schemes 

The previous section demonstrated how to view zero temperature ma- 
trix models as finite difference schemes for differential operators. Because 
the matrices were not random, the differential operators were not random 
either. This section extends to the finite (3 case, when randomness appears. 

6.1. Soft edge. 

6.1.1. Hermite — s- Airy. 

Definition 6.1. Let h = nT^/^. The n-by-n (3-Hermite matrix model 
scaled at the soft edge is 
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with DH^D ^ denoting the matrix in (|4.1() . 

The eigenvalues of -ff^ft display soft edge behavior as n ^ oo. The 
underlying reason, we claim, is that the matrix is a discretization of the 
stochastic Airy operator. The next theorem interprets iJ^fj as a finite 
difference scheme with mesh size h = n^^/"^ and grid points xi 
i = 1, . . . , n. 



— hi, 



Theorem 6.2. We have 



H 



soft 



H 



soft 



VP 



w, 



in which 



W 



Here, h 



1 

Vh 

-1/3 



Gi 

^(n-l)/3 



G2 

A(„-2)/3 



Gn 

, . . . , Gn are standard 

Gaussian random variables, and x^^^_j'^i3 ~ ^207i ^xfn-j)p ~ ("-^ j)/^)- -fwr- 



W has independent entries, Gi 



0{h ), uniformly 



ther, xfn-j)f3 ^'^^ mean zero and standard deviation 1 
for j such that Xj = hj < M , where M > is fixed. 

Proof. The derivation of the expression for -ff^f^ is straightforward. The 
mean of xfn-j)p exactly 0, and the variance is exactly 1 — h'^xj. □ 

We claim that the matrix W discretizes white noise on the mesh from 
Theorem 15. 21 The increment of Brownian motion over an interval {x,x + h] 
has mean and standard deviation ^/h, so a discretization of white noise 
over the same interval should have mean and standard devation . The 
noise in the matrix W has the appropriate mean and standard deviation. 

6.1.2. Laguerre Airy. 

The 2?T,-by-2?T, (3-Laguerre matrix model 



-1/3 



^soft 



DlPI 



Definition 6.3. Let h = (27i) 
scaled at the soft edge is 

{lp-'^Y 

in which Dl is the matrix D of Lemma [4.21 and Pl is the matrix P of the 
same lemma. The (2?! + l)-by-(27i + 1) (3-Laguerre matrix model scaled at 
the soft edge is 

2 / ^ r M^^- 



M: 



soft 



PmDiI - V2h-^/^I 
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in which D]\j is the matrix D of Lemma |4 . 31 and Pm is the matrix P of the 
same lemma. 

The eigenvalues of if^f^ and M^^^^ near zero display soft edge behavior as 
n oo. The underlying reason, we claim, is that the matrices themselves 
encode finite difference schemes for the stochastic Airy operator, as the 
next theorem shows. 

Theorem 6.4. The 2n-by-2n and {2n + l)-by-{2n + 1) (3-Laguerre matrix 
models scaled at the soft edge satisfy 



2 



and 



with 



soft soft 



VP 



2 



^ diag_;^(X(^^„)^,X(„_i)^, X(a+n-l)/3'X(„-2)/3> ■ ■ 



~2 -2 ~2 \ 

• ■ • ' X(q+2)/3: X(a+l)/3Ji 

diag_i(x^^, X(a+n-l)/3' X(n-l)/3i X(a+n-2);9i ■ • ■ 

~2 ~2 ~2 ~2 \ 

■ ■ ■ 1 X2/3; X(a+l)/3' Xpj Xapji 

and h — (2ri)^^/^. All 2n — 1 suhdiagonal entries of Wl and all 2n suhdi- 



agonal entries of Wm are independent, with Xr denoting a random variable 
with distribution Xr ^ ^2fin ^Xr ~ ^) ■ 

The entries of Wl and Wm have mean approximately and standard 
deviation approximately Therefore, we think of Wl and Wm as dis- 
cretizations of white noise on the mesh from Theorem 15.41 The situation 
is very similar to that in Theorem 16. 21 so we omit a formal statement. 

6.1.3. Overview of finite difference schemes for the stochastic Airy oper- 
ator. In light of Theorems 15.21 15.41 16.21 and 16.41 we make the follow- 
ing claim, -ff^ft, ^foft' ^^'^ -^soft discretize the stochastic Airy operator 
— —-^ + X + "^^'i fo'" finite and infinite fi. Because the low eigen- 
values of -ff^ft, -^foft' ^'^'i -^soft show soft edge behavior, it is natural to 
expect that the eigenvalues of show the same behavior. 

Conjecture 6.5. The kth least eigenvalue of the stochastic Airy operator 
follows the kth soft edge distribution, with the same value for (3. 

The conjecture now appears to be a theorem, due to a proof of Ramirez, 
Rider, and Virag |14j . 
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6.2. Hard edge. 

6.2.1. Laguerre —* Bessel. 

Definition 6.6. The n-hy-n j3-Laguerre matrix model scaled at the hard 
edge is 

in which h — — r—rr and F and fl are defined as in Section 1231 

The least singular values of L^'^.^ display hard edge behavior as n — > cx). 
We claim that this can be understood by viewing the matrix as a finite 
difference scheme for the stochastic Bessel operator with type (i) boundary 
conditions. The next theorem demonstrates this, using the same mesh seen 
in Theorem 15. 61 

Theorem 6.7. Let L^'^^^ he a matrix from the n-hy-n [3-Laguerre matrix 
model scaled at the hard edge. Adopting the notation of jy^.5| ) and Theorem 
15.61 and setting gi — ~ j3xi/h gi, we have 

(1) Lit, - e^-L3e-^-. 

(2) e'^^ = exp (-^ E^:;' x~"\3k^)) . 

(3) 5i,...,52n-i ire independent, and, for any e > 0, the random 
variables g\e/h'\^ ■ ■ ■ 7 52n-i have mean 0{\/h) and standard devia- 
tion 1 + 0{h), uniformly. 

The point of (3) is that the sequence -^gi, ■ ■ ■ , -^ff2n-i is a discretiza- 
tion of white noise. Hence, the expression for e''' in (2) is a discretization 

of ll}{Xi). 

Proof. Conclusions (1) and (2) are simply restatements of facts from (I4.5|) . 
For (3), the independence of gi, . . . ,g2n-i was already established in the 
context of H4.5|l . For the asymptotic mean and standard deviation, use 
the fact that log(xr/\A') has mean — ilog(r-/2) + ^Y{r/2) ~ 0{r^^) and 
standard deviation + 0{r^'^^^) as r oo. □ 

Theorem 15 .61 establishes that L^^^^ can be viewed as a finite difference 
scheme for the classical Bessel operator, 

Theorem 16.71 extends this connection to the /? < oo case. Noting the 
apparent convergence in distribution 

/ ^ 2«-i \ 
exp — -1= XI ^k^^'^ih^) i^ix), 

\ ^' k=[x-xo]/h J 
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^hard appears to be a finite difference scheme for tlie stocliastic Bessel 
operator with type (i) boundary conditions, 

Definition 6.8. The n-by-(n + 1) f3-Laguerre matrix model scaled at the 
hard edge is 

in which h — — F and are defined in Section [TS 

The small singular values of Mj^^'^^ display hard edge behavior as n ^ oo, 
because, we claim, the matrix is a finite difference scheme for the stochastic 
Bessel operator with type (ii) boundary conditions. The next theorem 
demonstrates this, using the same mesh seen in Theorem 15 .81 

Tiieorem 6.9. Let M^hard ^ matrix from the n-by-{n + 1) [3-Laguerre 
matrix model scaled at the hard edge. Adopting the notation of and 
Theorem \5.y\ and setting Cji — — (3xi j h gi , we have 

(1) Mttd-e^-MZ':,e~^-^. 

(2) ^o^p[-^j:t,x~'/'{~g,Vh)). 

(3) gi, ■ ■ ■ ,g2n are independent, and, for any e > 0, the random vari- 
ables g^^p{^, . . . ,g2n-i have mean 0(\/7i) and standard deviation 
1 + 0(h), uniformly. 

The point of (3) is that the sequence -^gi, ■ ■ ■ , -^Jl2n is a discretization 
of white noise. Hence, the expression for e"** in (2) is a discretization of 

ll){xi). 

Proof. Conclusions (1) and (2) are simply restatements of facts from (|4.6|) . 
For (3), the independence of 5i,...,52n was already established in the 
context of H4.6|l . For the asymptotic mean and standard deviation, use 
the asymptotics for chi-distributed random variables from the proof of the 
previous theorem. □ 

Hence, M^'^^ can be viewed as a finite difference scheme for the sto- 
chastic Bessel operator, 

<^d - e^-""AC-e-^=- ^I^J^r' - 

6.2.2. Jacobi — > Bessel. 

Definition 6.10. The n-by-n P-Jacobi matrix model scaled at the hard 
edge is 

jl3,a,b i_pD/3'a,6p 
•^hard ^^^22 ^' 
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in which h = 2n+a+b+i ^^"^ ^22° '' the bottom-right block of the 2n-by-2n 
/3-Jacobi matrix model jP^"-^^ . F and are defined in Section 

As n ^ 00, the small singular values of J^^^^ display hard edge behavior. 
We explain this fact by interpreting the rescaled matrix model as a finite 
difference scheme for the stochastic Bessel operator in Liouville normal form 
with type (i) boundary conditions. First, though, a lemma is required. 



Lemma 6.11. Suppose that 6 is a random angle in [0, ^] whose distribution 



is defined by cos^ ~ beta(c, d). Then log(tan^ 0) has mean ^^-^ — LIsl 



and variance r('j)r ^^J^^^ ^'^^ — 1_ r(c)r ('^^ , 

Proof, tan^ 6 = ^~™fg^ has a beta-prime distribution with parameters d, 
c. Hence, tan'^ 9 has the same distribution as a ratio of independent chi- 
square random variables, with 2c? degrees of freedom in the numerator and 
2c degrees of freedom in the denominator. Let X ~ X2d ^'^'^ ^ ^ xic be 
independent. Then the mean of log(tan^ 6*) equals i?[logX] — ^^[logF] = 
(log 2 -I- Yisy) — (log 2 -I- Yl^) — Yisy ~ ^"^^ tbe variance equals 

Var[logX] -t- Var[logy] = r(.)r'^(g-r'(.)- ^ r(c)r'V)-r'(c)- _ □ 

Theorem 6.12. Let J^^^'J^ be a matrix from the n-by-n (3-Jacobi matrix 
model scaled at the hard edge. Adopting the notation of \4-. 7| ) and Theorem 
\5.1U\ and setting 321-1 — — ■\/(/3x2j_i)/(2/i)i(log(tan^ 0i) — log(tan^^i)) 
and g2i = \J {f3x2i ) / (2/i) ^ (log(tan2 0'^) - log(tan2 9'^), we have 

/-.N jP.a.b D,,,„ jOO,a.b Dm 
\^) •-'hard ^ hard ^ 

(2) e''-' = e^pi^i^-^J^k'^i^ ^k^^^i9kVh)j in which R = -(log s^- 
logs„) if i ^ I, R ^ -(logs^_;^ - logs^_;^) - (logs„ - logs„) if 

i = 2j — \ is odd and greater than one, or R — (logSj — log Sj) — 
(logs„ — logs„) if i = 2j is even. 

(3) gi,...,g2n-i are independent, and, for any e > 0, the random 
variables g^^/h], ■ ■ ■ ,g2n-i have mean 0{\/}i) and standard devia- 
tion 1 -I- 0{h), uniformly. 

The point of (3) is that the sequence ■ • ■ j -^g2n~i is a discretiza- 

tion of white noise. Hence, the expression for e'^^ in (2) is a discretization of 
tpixi)^ . (The remainder term R has second moment 0{h) and is consid- 
ered negligible compared to the sum containing 2n~i terms of comparable 
magnitude.) 



Proof. Conclusion (1) is direct from (|4.7|l . 
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Now, we prove conclusion (2). According to (|4.7() . when i — 2j is even, 

n n—1 
d2j ^ ~ ^ (log Cfc - log Cfc) + ^ (log Sk - log Sk) 
k=j+l k=j 

n—1 n—1 

+ IZ(iog4 - logc;,) - ^(iog4 - logs'fc). 

k=j k=j 

Compare with 

-2n— 1 n n 



^ - ^/^{gkVh) ^ - ^ (log Cfc -log Cfc )+ ^ (logsfc-logsfc) 



n—1 n—1 



+ (log - 4) - II (log - log ^'fe^ 

The remainder term R is designed to cancel terms that occur in one ex- 
pression but not in the other. The argument for odd i is similar. 

The asymptotics in conclusion (3) can be derived from the explicit ex- 
pressions in the previous lemma. The details are omitted. □ 

6.2.3. Overview oj finite difference schemes for the stochastic Bessel oper- 
ator. Considering Theorems EHl EH [QHl O EH and WTR 

(1) i^ard discrctizcs with type (i) boundary conditions, for finite 
and infinite /3. 

(2) M^^^ discretizes with type (ii) boundary conditions, for finite 
and infinite /3. 

(3) Jhard'' discretizes jj^ with type (i) boundary conditions, for finite 
and infinite (3. 

Based on these observations and the fact that the small singular values of 
^hard' -^hard' ^nd J^ard^ approach hard edge distributions as n — + oo, we 
pose the following conjecture. 

Conjecture 6.13. Under type (i) boundary conditions, the kth least sin- 
gular value of the stochastic Bessel operator follows the kth hard edge dis- 
tribution with parameters (3, a. Under type (ii) boundary conditions, the 
hard edge distribution has parameters (3, a -\- \. This is true both for the 
original form, , and for Liouville normal form, . 



7. Numerical evidence 

7.1. Rayleigh-Ritz method applied to the stochastic Airy opera- 
tor. This section provides numerical support for the claim that stochastic 
Airy eigenvalues display soft edge behavior. Up until now, our arguments 
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have been based on the method of finite differences. In this section, we use 
the Rayleigh-Ritz method. 

To apply Rayleigh-Ritz, first construct an orthonormal basis for the 
space of L^((0,cx))) functions satisfying the boundary conditions for . 
The obvious choice is the sequence of eigenfunctions of . These func- 
tions are Vi{x) — Ai(— x -I- Q), i — 1,2,3,..., in which is the 
ith zero of Airy's function Ai. Expanding a function v in this basis, 
V = civi + C2V2 + ■ ■ ■ , the quadratic form (v^A'^v) becomes 

{V,A'^V)^ V ({c,V„A°°CjVj) + ^ f (c,Vi){CjVj)dB 

,^j->i V VP Jo 



Note that the stochastic integral viVjdB is well defined, and its value 
does not depend on specifying an Ito or Stratonovich interpretation, be- 
cause Vi and Vj are well behaved and not random. (In fact, the joint 
distribution of the stochastic integrals is a multivariate Gaussian, whose 
covariance matrix can be expressed in terms of Riemann integrals involv- 
ing Airy eigenfunctions.) Introducing the countably infinite symmetric K, 

2 

7^. 



K = { -CAj + -j= / ViVjdB 



i,i=l,2,3,... 

the quadratic form becomes c^Kc, in which c = (ci, C2, C3, . . . )^, the vector 
of coefRcients in the basis expansion. 

According to the variational principle, the least eigenvalue of A^ equals 
mi\\y\\^i{v,A^v), which equals min||j,||=i c^i^c, which equals the minimum 
eigenvalue of K. This suggests a numerical procedure. Truncate K, taking 
the top-left l-hy-l principal submatrix, and evaluate the entries numerically. 
Then compute the least eigenvalue of this truncated matrix. This is the 
Rayleigh-Ritz method. 

The histograms in Figure ^21 were produced by running this procedure 
over 10^ random samples, discretizing the interval (0,86.9) with a uniform 
mesh of size 0.05 and truncating K after the first 150 rows and columns. 
The histograms match the soft edge densities well, supporting the claim 
that the least eigenvalue of A^ exhibits soft edge behavior. 

7.2. Rayleigh-Ritz method applied to the stochastic Bessel op- 
erator. Now consider applying the Rayleigh-Ritz method to the stochasic 
Bessel operator in Liouville normal form with type (i) boundary conditions. 
Liouville form is well suited to numerical computation because the singular 
functions are well behaved near the origin for all a. We omit consideration 
of type (ii) boundary conditions for brevity. 
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Two orthonormal bases play important roles, one consisting of right sin- 
gular functions and the other consisting of left singular functions of J"^ , 
from H3.9|l . For i = 1, 2, 3, . . . , let Vi{x) be the function \fxia{t,ix)^ normal- 
ized to unit length, and let Ui(x) be the function ■\/a;ja+i(^ia;), normalized 
to unit length, in which is the ith zero of ja- 

The smallest singular value of jj^ is the minimum value for ^^"1"° | "^^ = 

^— — '"^y^y' Expressing v bs v ~ ftp^ and expanding / in the basis 

vi,V2,V3, ... as f — X^i^i '^j'^^i' ttic norm of ||j7^w|| becomes 

.1/2 




In terms of the countably infinite symmetric matrix 



we have || J7jf = {c^ Kc)^/'^ , in which c = (ci, C2, C3, . . . . For the norm 
of w, we find 



11^1 = 11/^^11 = 




in which c is defined as above and M is the countably infinite symmetric 
matrix 



M ^ { I ^'^^ViVjdt 



The least singular value of equals min , which equals min {jtt^j^ 

which equals the square root of the minimum solution A to the generalized 
eigenvalue problem Kc = XMc. To turn this into a numerical method, 
simply truncate the matrices K and M, and solve the resulting generalized 
eigenvalue problem. 

The histograms in Figure [T31 were produced using this method on 10^ 
random samples of the stochastic Bessel operator, discretizing the interval 
(0, 1) with a uniform mesh of size 0.001 and truncating the matrices K and 
M after the first 75 rows and columns. The histograms match the hard 
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edge densities well, supporting the claim that the least singular value of 
the stochastic Bessel operator follows a hard edge distribution. 

7.3. Smoothness of eigenfunctions and singular functions. Up to 

this point, we have proceeded from random matrices to stochastic oper- 
ators. In this section, we reverse direction, using stochastic operators to 
reveal new facts about random matrices. Specifically, we make predictions 
regarding the "smoothness" of Hcrmite eigenvectors and Jacobi CS vectors, 
using the stochastic operator approach. Verifying the predictions numeri- 
cally provides further evidence for the connection between classical random 
matrix models and the stochastic Airy and Bessel operators. 

First, consider the eigenfunctions of the stochastic Airy operator. The 
fcth eigenfunction is of the form in which fk € C^((0,oo)) and G 
C^/^~((0,oo)) is defined by (|3.2|1 . In hght of the claim that H^^^^ encodes 
a finite difference scheme for , the kth eigenvector of H^^^^ should show 
structure indicative of the fcth eigenfunction of . For a quick check, 
consider the ratio of two eigenfunctions/eigenvectors. The kth eigenfunc- 
tion of A^ is of the form fk4>, which does not have a second derivative (with 
probability one) because of the irregularity of Brownian motion. However, 
the ratio of the fcth and ^th eigenfunctions is = which, modulo poles, 
has a continuous second derivative. Therefore, we expect the entrywise ra- 
tio between two eigenvectors of i?f,ft , L^^it ' -^soft "smoother" than 
a single eigenvector. Compare Figure [T^ 

Next, consider the singular functions of the stochastic Bessel operator in 
Liouville normal form. The fcth right singular function is of the form 
and the fcth left singular function is of the form gkip"^, in which fk,gk G 
Ci((0, 1)) and ?A G Ci/2-((o, 1)) is defined by (gTI). The situation is similar 
to the Airy case. Any one singular function may not be differentiated in the 
classical sense, because of the irregularity of Brownian motion. However, 
the ratio of two singular functions is smooth. We expect the singular vectors 
of -^hard' -^hard' ^havd show similar bchavior. Compare Figure [T31 

8. Preview of the stochastic sine operator 

We have seen that the eigenvalues of the stochastic Airy operator display 
soft edge behavior, and the singular values of the stochastic Bessel operator 
display hard edge behavior. Is there a stochastic differential operator whose 
eigenvalues display bulk behavior? Because of the role of the sine kernel in 
the bulk spacing distributions, it may be natural to look for a stochastic 
sine operator. In fact, jl6| provides evidence that an operator of the form 



(8.1) 





_d_ ■ 

dx 


+ 


"noise" 


"noise" 


_d_ 

dx 






"noise" 


"noise" 
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may be the desired stochastic sine operator. This operator is discovered 
by scahng the Jacobi matrix model at the center of its spectrum, and an 
equivalent operator, up to a change of variables, is discovered by scaling 
the Hermite matrix model at the center of its spectrum. The exact nature 
of the noise terms in (|8.1|l is not completely understood at this point. A 
change of variables analogous to those that transform to (|3.3|l and 

1)3.6(1 to 1)3. 8|l would be desirable. 

Acknowledgement. The authors thank loana Dumitriu for the inspiration 
provided by her ideas. Figure [TTTl was produced with the help of Per-Olof 
Persson's software. 



References 

loana Dumitriu and Alan Edelman. Matrix models for beta ensembles. J. Math. 
Phys., 43(ll):5830-5847, 2002. 

loana Dumitriu and Alan Edelman. Eigenvalues of Hermite and Laguerre ensembles: 
large beta asymptotics. Ann. Inst. H. Poincare Probab. Statist., 41(6):1083— 1099, 
2005. 

Alan Edelman. SIAM Conference on Applied Linear Algebra. The College of 
William and Mary, Williamsburg, VA, July 2003. 

Alan Edelman and Brian D. Sutton. The beta-Jacobi matrix model, the CS decom- 
position, and generalized singular value problems. Preprint. 

P. J. Forrester. The spectrum edge of random matrix ensembles. Nuclear Phys. B, 
402(3):709-728, 1993. 

P. J. Forrester. Painleve transcendent evaluation of the scaled distribution of the 
smallest eigenvalue in the laguerre orthogonal and symplectic ensembles. 2000. 
http : //arxiv . org/abs/nlin . SI/0005064. 

Michio Jimbo, Tetsuji Miwa, Yasuko Mori, and Mikio Sato. Density matrix of an 
impenetrable Bose gas and the fifth Painleve transcendent. Phys. D, 1{1):80— 158, 
1980. 

Kurt Johansson. Shape fluctuations and random matrices. Comm. Math. Phys., 
209{2):437-476, 2000. 

Iain M. Johnstone. On the distribution of the largest eigenvalue in principal com- 
ponents analysis. Ann. Statist., 29(2):295-327, 2001. 

Rowan Killip and Irina Nenciu. Matrix models for circular ensembles. Int. Math. 
Res. Not., (50):2665-2701, 2004. 

A. B. J. Kuijlaars and M. Vanlessen. Universality for eigenvalue correlations from 
the modified Jacobi unitary ensemble. Int. Math. Res. Not., (30):1575-1600, 2002. 
Madan Lai Mehta. Random matrices. Academic Press Inc., Boston, MA, second 
edition, 1991. 

Bernt 0ksendal. Stochastic Differential Equations: An Introduction with Applica- 
tions. Springer- Verlag, Berlin, fifth edition, 2000. 

Jose A. Ramirez, B. Rider, and Balint Virag. Beta ensembles, stochastic airy spec- 
trum, and a diffusion. Preprint. 

Jack W. Silverstein. The smallest eigenvalue of a large-dimensional Wishart matrix. 
Ann. Probab., 13(4):1364-1368, 1985. 

Brian D. Sutton. The stochastic operator approach to random matrix theory. PhD 
thesis, Massachusetts Institute of Technology, Cambridge, MA 02139, June 2005. 



FROM RANDOM MATRICES TO STOCHASTIC OPERATORS 



41 



[17] Gabor Szego. Orthogonal polynomials. American Mathematical Society, Providence, 
R.I., fourth edition, 1975. American Mathematical Society, Colloquium Publica- 
tions, Vol. XXIII. 

[18] Craig A. Tracy and Harold Widom. Introduction to random matrices. In Geometric 
and quantum aspects of integrable systems (Scheveningen, 1992), volume 424 of 
Lecture Notes in Phys., pages 103—130. Springer, Berlin, 1993. 

[19] Craig A. Tracy and Harold Widom. Level-spacing distributions and the Airy kernel. 
Phys. Lett. B, 305:115-118, 1993. 

[20] Craig A. Tracy and Harold Widom. Level-spacing distributions and the Airy kernel. 
Comm. Math. Phys., 159(1);151-174, 1994. 

[21] Craig A. Tracy and Harold Widom. Level spacing distributions and the Bessel 
kernel. Comm. Math. Phys., 161(2):289-309, 1994. 

[22] Craig A. Tracy and Harold Widom. On orthogonal and symplectic matrix ensem- 
bles. Comm. Math. Phys., 177(3):727-754, 1996. 

[23] Craig A. Tracy and Harold Widom. The distribution of the largest eigenvalue in the 
Gaussian ensembles: /3 = 1,2,4. In Calogero-Moser-Sutherland models (Montreal, 
QC, 1997), CRM Ser. Math. Phys., pages 461-472. Springer, New York, 2000. 

[24] Hale F. Trotter. Eigenvalue distributions of large Hermitian matrices; Wigner's 
semicircle law and a theorem of Kac, Murdock, and Szego. Adv. in Math., 54(1):67- 
82, 1984. 

Alan Edelman, Department of Mathematics, Massachusetts Institute op 
Technology, Cambridge, MA 02139 

Brian D. Sutton, Department of Mathematics, Randolph- Macon College, 
Ashland, VA 23005 



